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ABSTRACT 

An earlier paper gave solutions for the mean time rates of change of 

i 

orbital elements of satellite atoms in an exosphere influenced by solar 
radiation pressure. Each element was assumed to behave independently. Here 
the instantaneous rates of change for three elements (e, ft, and = 01 + n) are 
integrated simultaneously for the case of the inclination i = 0 . The results 
(a) confirm the validity of using mean rates when the orbits are tightly bound 
to the planet and (b) serve as examples to be reproduced by the complicated 
numerical solutions required for arbitrary inclination. Strongly bound 
hydrogen atoms escaping from Earth due to radiation pressure do not seem a 
likely cause of the geotail extending in the anti-sun direction. Instead, 
radiation pressure will cause those particles’ orbits to deteriorate into the 
- Earth's atmosphere. Whether loosely bound H atoms are plentiful enough to 
create the geotail depends on their source function ve rsus r; that question is 
beyond the scope of this paper. 
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INTRODUCTION 

In an earlier paper (Chamberlain, 1979, hereinafter referred to as "Paper 
I") we developed equations for the mean rates of change of the orbital 
elements of satellite particles In a planetary exosphere that is subjected to 
solar radiation pressure. The results showed a surprising tendency for the 
perigees of direct orbits to lock into stable positions westward of the planet 
as seen from the sun (and eastward for retrograde orbits). These stable 
positions of perigee are also close to the positions where perigee is most 
rapidly lowered, vacating the orbits (Paper I, Fig. 4). 

It would thus appear that satellite orbits should be efficiently depleted 
in time scales of the order of [Paper I, Eq. (17)] 

T 0 = (y/a) 1/2 /f 
= 1 J055 x 10 6 /(a/R E ) 1/2 sec, 


where y = GM^, f is the acceleration due to radiation pressure, and the 
orbital notation is conventional and follows that of Paper I. To pursue this 
matter further, we have investigated the simultaneous solution of the elements 
<{> s (longitude of perigee from the sun), e (eccentricity) and (longitude 
from the sun of the ascending node). In these solutions we set the inclina- 
tion i = 0. For the last two elements we will use mean values. The mean rate 
of change of the semi -major axis vanishes: <da/dt> = 0. And the time of 
perigee passage, x, offers special problems that make it desirable to adopt 
its mean value. 
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INSTANTANEOUS AND MEAN RATE EQUATIONS 

The equations for perturbed elements are given in terms of the external 

force components by Burns (1976), Burns et al. (1979), and in Paper I. In 

I 

Paper I, by the way, the value of de/dt is incomplete, since it was derived by 
holding a. - const., in anticipation of averaging over the orbit. In reality a. 
varies over s perturbed orbit (only its mean time derivative vanishes) and the 
correct Eq. (10) of Paper I is 

de/dt = [a(l - e 2 )/u] 1/2 [ - f T (l - e 2 )/e(l + 3 cos v) 

+ fp sin v + fy(l + e cos v)/e]. (2) 

The error does not affect the derived mean value, <de/dt>, or any other 
- results of Paper I. 

To integrate these rate equations, the force components must be expressed 
in terms of the particle's position in space [Paper I, Eqs. (2), (6), and 
(7)]. The instantaneous equations are then, for the semi-major axis, 

da/dt = 2f [a 3 /y(l - e 2 )] 1/2 [cos sin (u> + v) 

s 

+ sin cos (to + v) cos i + e sin n cos to cos i 

+ e cos sin to], (3) 


for the eccentricity, 
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de/dt * f[a(l - e 2 )/u] 1/2 { - [(1 - e 2 )/e(l + e cos »)] [cos fi g sin (<d + v) 

+ sin cos (id + v) cos i] - sin v [cos fi g cos (w + v) 

- sin fi $ sin (id + v) cos ij + [(1 + e cos v)/e] x 

[ cos fig sin(w + v) + sin fi $ cos (<d + v) cos i]}, (4) 

for the longitude of the ascending node, 

dfi/dt = - f [a(l - e 2 )/y] 1 ^ 2 sin fi sin (» + v)/(l + e cos v), (5) 

S *5 

for the inclination, 


di/dt = - f [a (1 !- e 2 )/y]^ 2 sin fi $ sin i cos (id + v)/ 

(1 + e cos v), 


( 6 ) 


for the argument of perigee, 

dai/dt + cos i dfig/dt = f [a(l - e 2 )/ye 2 ]^“ x 

[cos fi cos (<d + v) cos v - sin fi sin (<d + v) cos i cos v 
s s 

* 

+ cos fig sin (id + v) sin v (2 + e cos v)/(l + e cos v) 

+ sin fig cos (<d + v) cos i sin v (2 + e cos v)/(l + e cos v) ] , (7) 



and finally (Burns, 1976) for the rate of change of the time of perigee 
passage, 

dt/dt » f{3(r - t)[a/u(l - e^)]*^ e sin v - [a^(l - e^)/ue][cos v - 

2e/(l + e cos v)]}[ - cos ft g cos (w + v) + sin « s sin (u> + v) cos i] 

+ f {3(t - t)[a/u(l - e^)]^(l + e cos v) 

+ [a <: (l - e c )/ue][sin v(2 + e cos v)/(l + e cos v) ] } x 

[cos ft g sin (ti> + v) + sin ft g cos ( u> + v) cos i ]. (8) 

-In these equations v is the true anomaly and ft g is related to ft (the longitude 
of the node measured -from the vernal equinox) by ft = ft - X, where 
X = n $ (t *■ t ) and n $ is the mean solar motion about the planet. 

In all these equations we need v(t) to carry out the integrations. This 

functional relationship comes from Kepler's equation, 

M = n p (t - t ) - e - e sin e (9) 

where M is the mean anomaly, e the eccentric anomaly, n p the mean motion of 
the particle, and where 
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Kepler's equation may be solved by Newton-Raphson iteration to find 
v(t) when t(t) is known, x(t) is found in turn from Eq. (8) if v(t) is 
known. Thus (8), (9), and (10) have to be solved as a simultaneous set. 

An additional complication is that t - t itself appears in the integrand 

* 

f 2 

of Eq. (8). From the law of areas, dv/ot = H/r , we have 

t(v) - t - [a 3 (l - e 2 ) 3 /u] 1 / 2 0 J y dv/(l + e cos v ) 2 

= [a 3 ( 1 - e 2 ) 3 /u]^ 2 { - e sin v/(l - e 2 )(l + e cos v) 

+ 2(1 - e 2 )" 3/2 tan _1 [(l - e) 1/2 (1 + e)" 1/2 tan v/2]} (11) 

To avoid the simultaneous solution of Eqs. (8) [with (11)3, (9), and 
(10), we average Eq. (8) over a cycle and then use <x(t)> in lieu of t(t) in 
_ solving Kepler's equation. For i = 0 the force terms in the rate equations 
[Paper I, Eqs. (2), (6), and (7)] simplify to sin (<j> s + v) and cos ($ s + v), 
where 4> s is the longitude of perigee from the sun. The mean rate for x is 

<dx/dt> =(3a 2 f cos <j> s /2ye)[l + 2e 2 

+ 2(1 ; g2) e 0 J 2 * F(v)dv], (12) 

where 


p( v ) - sin v tan~* {[(1 - e) / ( 1 + e)]*^ tan v/2> 
(1 + e cos v) 2 


( 13 ) 
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For small e, Eq. (12) is comparable to 

<d* s /dt>/n p * (3a 2 f cos 4> s /2ye)(l - e 2 ) 1/2 , (14) 

t, 

» 

where n p is the mean rate of motion of a particle about the plnnet 
(n p * 2-n/P * /a 3 ^ 2 ). That Eqs. (12) and (14) must agree in the limit of 

e 0 follows from the consideration that the actual angular rate of particle 
motion, dv/dt, is equivalent to a^/At and to the mean particle rate, n p , as 
the orbit approaches circularity. 

The integral of Eq. (13) may be evaluated by Gaussian integration and 
then represented as a power series in e. This procedure gives 

<dr/dt> = (3 a 2 f cos <J> s /2ye)[l + 2e 2 

+ |(1 - e 2 ) e(0.968 + 7.30e - 40.482 e 2 

+ 79.0617 e 3 - 43.6806 e 4 )]. (15) 


Figure 1 illustrates the accuracy of this polynomial representation of F(v). 

The integrals to obtain the mean rates may be found from the tables of 
Gradshteyn and Ryzhik (1965, Sect. 2.55). They are of the form 


and the identity 


(n) ,2n n 

I = / dv g(v)/(l + e cos v) , 

g o 


If/ 0 " ^ " e ^ cos v = ij 

g g g 


(16) 

(17) 


is occasionally useful. 
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SOLUTION AND CONCLUSION 

With i * 0 and <da/dt> * 0 and with <x(t)> used to solve Kepler's 
equation for v(t), we can readily integrate Eqs. (5), (6), and (7) simultan- 
eously. The results are shown in Figure 2 along with solutions for the mean 
rates obtained from the equations of Paper I. 

The similarity between mean and instantaneous solutions confirms the 
validity of using the mean rates of change when the orbits are tightly bound 
to the planet. In paper I [Eq. (18)] the criterion for using means rates was 
specified as 


^ = 2.086 x 10 2 /(a/R E ) 2 » 1. (18) 

In the present case, Eq. (1) gives T 0 = 4.7 x 10 5 sec = 5.5 days, and we 
have T Q /p = 8.3, so that the criterion is crudely satisfied. 

The eccentricity for a particle collision with the planet is 

e CQll = 1 - l/(a/R E ). (19) 

From Figure 2 this condition is met at e co i -j = 0.8 at the end of 5 orbits or 
about 3.3 days, which (as expected) is the order of T Q . The decay of eccen- 
tricity shown in Figure 2 is a clear illustration of the orbital instability 
first noted in Paper I: The longitude of perigee (for a direct orbit) moves 

asymptotically towards 90°, which is the region where perigee is most rapidly 
lowered (see Fig. 3). The combined effect is to rapidly vacate orbits in the 
ecliptic plane (or in any other orbit whose plane contains the Earth- Sun 
line). 
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Figures 3, 4, 5, 6, and 7 show sample calculations for loosely bound 
orbits. In the extreme case of a/R^ » 1, radiation pressure can remove bound 
orbits within a fraction of an orbit, but in intermediate cases (a/R^ 10) the 
atoms may be either removed or forced into collision with the planet as the 
eccentricity is increased. Further, the escaping orbits show no overwhelming 
preference for exit in the anti-sun direction. The existence of a geotail, 
observed in Ly a, may be due to charge exchange of loosely bound or escaping 
atoms with magnetospheric, high-energy protons, followed by a radiation- 
pressure impulse. However, securely bound atoms ( i . e . , ones that are only 
slightly perturbed during one revolution) could not be perturbed sufficiently 
to cause the observed anti-solar asymmetry to the geocorona. 
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FIGURES 

Fig. 1. Polynomial representation of the integral of Eq. (13). The solid 
curve is 

[2(1 - e 2 )e/n) / 2 V(v)dv 
o 

and the dashed curve gives [4(1 - e 2 )e/*]( 0.968 + 7.30 e 

- 40.482e 2 + 79.0617e 3 - 43.6806 e 4 ). 

Fig. 2. Comparison of the instantaneous elements, e, $ s ( * w), and 

(dotted curves) with their mean values averaged over an orbit (solid 
lines) for satellite H atoms in Earth orbit. The initial elements are 
a = BR^, i * 0, e * 0.1, 4> s a 0, and & s * 45°. 

Fig. 3. Orbits of an H atom at a/R e = 5. Note that the perigee, initially in 
the sunward direction, moves to 4> s * 90° by the fifth orbit. On the 
sixth orbit, where e > 0.8, the particles will crash into the planet. 

Flo, 4. Progression of the eccentricity for orbits at a = 10 and 15 Earth- 
radii. The initial values are e Q = 0.5. These values are actual 
eccentricities for continuous radiation pressure; the mean elements would 
be a poor approximation. Quantum accelerations, at mean thrust intervals 
of 435 sec., produce some dispersion about the continuum acceleration. 

Fig. 5. Progression of the perigee longitude from the sun, $ , for orbits at 
a =10 and 15 Earth-radii. See legend to Fig. 4. 
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Fig. 6. Orbit of an H atom about Earth with a/P. fi * 10. obtained from the 

perturbed elements. The progression of e and $ s are given In Figs. 4 and 

5. This particle collides with the planet on its second orbit. J? it 
missed collision it would escape, but not in the anti-sun direction and 
not contributing to a geotail of luminescent H, 

Fig. 7. Orbit of an H atom about Earth with a/R e * 15, obtained from the 

perturbed elements. The progression of e and $ are given in Figs. 4 and 

5. This particle escapes on its first orbit, in a direction that does 
not contribute to a geotail of luminescent H. 
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